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ABSTRACT 

We compute the skewness t% and the corresponding hierarchical amplitude T3 
of the divergence of the velocity field for arbitrary non-Gaussian initial conditions. 
We find that T3 qualitatively resembles the corresponding hierarchical amplitude for 
the density field, S3, in that it contains a term proportional to the initial skewness, 
which decays inversely as the linear growth factor, plus a constant term which differs 
from the corresponding Gaussian term by a complex function of the initial three- and 
four- point functions. We extend the results for S3 and T3 with non-Gaussian initial 
conditions to evolved fields smoothed with a spherical tophat window function. We 
show that certain linear combinations, namely S3 + 5T3, 53 + T3, and S3 + £3, lead to 
expressions which are much simpler, for non-Gaussian initial conditions, than S3 and 
T3 (or S3 and £3) considered separately. 



Subject headings: cosmology: theory - galaxies: clustering - large-scale structure of 
Universe 



Introduction 



The standard cosmogenesis lore attributes the formation of large-scale structure to the 
enhancement of primordial density fluctuations by gravity. Given that observations on scales 
larger than 10 h _1 Mpc show the amplitude of the rms fluctuations to be less than unity, one can 
successfully apply perturbative techniques to follow the evolution of the initial probability density 
function (hereafter PDF). However, a complete analysis of the problem requires some knowledge 
of the statistical nature of the initial fluctuations. Simple inflationary scenarios naturally produce 
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initial density fluctuations characterized by a Gaussian PDF, whereas models based on cosmic 
strings, global texture, or inflation models with multiple scalar fields (Vilenkin 1985; Kofman 
1991; Gooding et al. 1992) provide fluctuations which are initially non-Gaussian. 

A great deal of analytic work has been done for Gaussian initial conditions. In particular, the 
hierarchical amplitudes corresponding to the skewness S3 and the kurtosis S4 have been known 
for some time (Peebles 1980; Fry 1984), and more recently Bernardeau (1992) has provided a 
formalism for the calculation of the full hierarchy. This work was extended to the skewness of 
the velocity field by Bernardeau (1994a) and Bernardeau et al. (1995), and some progress has 
been made toward estimating the evolution of the full PDF (Kofman et al. 1994; Bernardeau & 
Kofman 1995; Protogeros k, Scherrer 1996). However, a firm conclusion has not yet been reached 
as to the nature of initial conditions that would generate, after evolution, large scale structures 
with the same statistical characteristics as the ones obtained from observations such as the QDOT 
survey or the 1.2 Jy IRAS survey (Coles & Frenk 1991; Bouchet et al. 1993). One reason is that 
although the statistical description of the evolved fields has been calculated in detail for Gaussian 
initial conditions, as noted above, the same cannot be said for non-Gaussian fields. 

For the case of non-Gaussian initial conditions, a general expression for the evolution of 
the skewness was derived by Fry & Scherrer (1994), and for the kurtosis by Chodorowski & 
Bouchet (1996). Here we extend this earlier work by examining the evolution of the divergence 
of the velocity field {6 = V • v/aH) for general non-Gaussian initial conditions. Specifically, we 



denote statistical indicators related to the density field and T to denote the ones related to the 
divergence of the velocity field.) As noted above, the skewness of 6 for Gaussian initial conditions 
has already been calculated (Bernardeau 1994a; Bernardeau et al. 1995; Bernardeau et al. 1996), 
as well as the correlation between 5 and 6 in the quasi- linear regime (Chodorowski & Lokas 1996). 

The calculation of T3 and £3 is presented in the next section. As in the case of the skewness 
of the evolved density field, we find three contributions to T3: i) the initial, linearly-evolved 
skewness, which decays as the inverse of the linear perturbation growth rate, ii) a "Gaussian" 
contribution, identical to the value of T3 for Gaussian initial conditions, and constant in time, and 
iii) a term, also constant in time, which is a complex function of the initial three- and four-point 
correlations. The values of T3 and S3 for the smoothed final density field are calculated in Section 
3. We discuss our results in Section 4 and show that certain linear combinations of T3 and S3 give 
simpler expressions than each considered separately. 



calculate the normalized skewness £3 = — — /t-, as well as the corresponding hierarchical coefficient 
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T3 = 2 for arbitrary non-Gaussian initial conditions. (Throughout the paper, we use S to 





2. Calculation of T 3 
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In what follows we adopt Peebles's (1980) notation in a line of argument paralleling that of 
Fry and Scherrer (1994). To simplify our derivation, we define 9 = H9 = V ■ v/a. We begin by 
expanding 5, v, and 9 perturbatively in the form: 



5 = 5^+6^ + ... 

v = v (l) +v (2) + ... 

9 = e (1) + 9 (2) + ... 



so that the skewness of 9 is given by 

6(0) = (o 3 



+ 3 



(i) 



(2) 



The first-order and second-order solutions for 5 are the well-known linear result: 

5^ = D 1 (t)5(x,t ), (3) 
where D\{t) is the usual growing mode solution, and (Peebles 1980), 

5 {2) = pi - 5 0yi A 0A + |Ao,ijAo,ij, (4) 

where we have defined 5q = 6^ , and Ao is defined to be: 

A„ = i/ rf V«x<)^. (5) 

To go from these expressions to the values for 6^ = V • v^ 2 ) /a and 9^ = V • v^/a, we 
substitute the expansions for 5 and v into the continuity equation: 

j + i V • (1 + <5)v = 0, (6) 
leading to the linear and second-order equations: 

jW + Iv-yW = 0, (7) 
a 

( 5(2) + I v . v ( 2 ) + i<5( 1 )V-v( 1 ) + iv( 1 ) •VJW = 0. (8) 
a a a 

The first-order equation gives nf* = at)\/S.Q^ and V • = — aZ>i<5o(x), so that 

= -DkJo(x), (9) 

while the second-order equation gives: 

0(2) = _ Dl JJ l (1 6 2 _ 6oiAqi + * Ao . jAo .^ . (10) 
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From equations (B) and fli~0|), we obtain: 



and 



6f(D V §(2) 



0(D 

3 
7 



-DH[6< 



(11) 



(12) 



We now simplify the various terms in equation (|12| ) using the results of Peebles (1980) and Fry & 
Scherrer (1994). 

The second-, third-, and fourth-order moments of the density field can be expressed in terms 
of the irreducible 2-, 3- and 4-point correlation functions £12, Ci23 5 an d 7/1234: 



<<5(xi)<5(x 2 )) 
(ff(x 1 )5(x 2 )<y(x 3 )) 
(<5(x 1 )5(x 2 )5(x 3 ) ( 5(x4)) 



62, 

Cl23; 

^12^34 + Cl3^24 + ^23^14 + ??1234- 



(13) 



where all the moments are functions of |x — x'| due to the assumed homogeneity and isotropy of 
the density field distribution. Therefore, 

:<5o) = 3£o 2 (0)+r/o(0), (14) 



and 



1 

4tt 



d 3 x'- : 



1 



+ 5S ' rfV 



X — X 

1 



X — X 



7i,i£o(0)&(x-x'),i 
-,i?7o(x,x,x,x'),i, 



(15) 



where the zero subscripts indicate linearly-evolved quantities. The first term in equation (15) 
yields £,o(0), whereas the second term yields §?7o(0) upon integration by parts, using the fact that 
r/(x, x, x, x) = 7/(0) from isotropy. Thus, 



<^o,iAo,i) =£g(0) + -770(0). 



Finally, we evaluate the last term in equation (|12|): 



1 



d 3 x'd 3 x"- 



1 



1 



x — x' x — x" 



<5 2 (x)<5 (x')<5o(x")). 



(16) 



(17) 



(4tt) 2 

The fourth moment on the right-hand side can be expressed as 

;5 2 (x),5o(x / )<5o(x")) = £o(0)£o(x' - x") + 2£ (x - x%(x - x") + t/ (x, x, x', x"). (18) 



The first two terms in the double integration in equation (|i~7|) yield £q(0) an( ^ §Co (0) respectively, 
for a total contribution of §£o(0)- The last term yields I[r)o] + |r/o(0), where we define (Fry & 
Scherrer 1994): 



(47T) 



dVdV'r/o(0,0,x',x 



, „,6P 2 (x'-x") 



(19) 
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with P2(x' ■ x") the second Legendre polynomial. Therefore, 



<5gAo.ij-Ao.ij) = -Co(0) + 3%(0). 



(20) 



Combining the results of equations (|I4|), (|l6|), and (20), equation (12) becomes 



§<n 2 #2) 



D\Di 



^2(0) + ^-ijb(0) + 



(21) 



Substituting equations (|ll]) and (pi|) into equation @, we obtain our final expression for Cg(0) : 



C e -(o) = -b\D x 



^ + feo 2 (o) + ^o(o ) + ^/N 



(22) 



To calculate the hierarchical amplitude T3 or the normalized skewness £3, we must also 
calculate £g(0) = /> which can be derived in a calculation similar to that for Q. We have 



with 



and 



0(1) 



^ 2 e (o) 



gf(i)gl(a)) = DfA < -<5 3 - 5 5o, 4 A , 4 + -^Aq^Ao,^ 



' 3 4 
-^o - 5 5o,iAo,i + - 

In a manner similar to our previous derivations, we obtain 



[4 

(^oAo.ijAo.ij 



Co(0), 

§Cb(o), 

/[Co] + ~Co(0), 



(23) 

(24) 
(25) 

(26) 
(27) 
(28) 



where /[Co] is the integral analogous to equation (|i~9|), integrated over Co(0,x',x"), rather than 



'/o 



(0,0,s',x"). Therefore, {9^9^ 



1 



-Co(0) + -I[Co] 



and 



b\D x 



(29) 



o.(-i 



Combining equations ( |2"2| ) and (|29|), with (9 = — 9, and defining /(f2) = 7777^ ~ 0° ". we 

— — ^ ^ 

obtain our expressions T3 and £3 for arbitrary non-Gaussian initial conditions. For the hierarchical 

amplitude T3, we get 

1 



m) 



1 Co(0) , 26 6^(0) 12 /M 



10C 2 (0) 16Co(0)/[Co] 



^(t)Co 2 (0) 7 7^(0) 7C 2 (0) 21^(0) 7 $(0) 



(30) 
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where only terms of order up to 0(a°) have been retained. In the case of Gaussian initial 
conditions, all of the terms except the second term are zero, yielding the standard result 
(Bernardeau 1994a, Bernardeau et al. 1995) 



T 3 G 



20 



-0.6 



(31) 



For comparison the hierarchical coefficient for the density field is (Pry & Scherrer 1994) 

1 Co(0) 34 10 7^(0) 6 I[r) ] 26 C 2 (0) 8 Co(0)/[Co] 
3 A£ 2 (0) 7 7£o 2 (0) 7£o 2 (0) 21£ 3 (0) 7 £g(0) " 



(32) 



We see that T 3 , like £3, contains three distinct contributions: a term proportional to the initial 
skewness which decays away as 1/D\{t), a "Gaussian" piece which is constant and identical to the 
hierarchical amplitude in the Gaussian case, and a third contribution, also constant in time, which 
is a complex function of the initial skewness, the initial kurtosis, and various integrals over the 
initial three- and four-point functions. In fact, a comparison of equations (32) and (|3C|) indicates 
that the terms in the two equations are identical functions of the initial density field; only the 
coefficients multiplying the various terms are different. We will exploit this fact in Section 4. 

For the normalized skewness t% we obtain 



h = -^3,0 

26 6 7/o(0) 12 JM _ .5_Co(0) _ 12 Co(0)J[Co 
7 + 7£ 2 (0) + 7£ 2 (0) 14^(0) 7 $(0) 



0"o(*), 



(33) 



which can be compared with the corresponding normalized skewness for the density (Fry & 
Scherrer 1994) 



S3 



+ 



S3,0 

34 

y 



10770(0) + 6/[r/ ] 



13C 2 (0) 



7 £o 2 (0) 7^(0) 14 £3(0) 



6 Co(0) J[Co] 

7 £ 3 (0) 



a (t). 



(34) 



In these expressions, ao(t) is the linearly-evolved rms fluctuation: crg(t) = £o(0), and 53^ is the 
(constant) linearly-evolved skewness: s^q = Co(0)/cg(t). Note that ^3, unlike T 3 , is independent of 



3. Smoothed T 3 and S3 results 



Although the results derived in the previous section are interesting from a formal point of 
view, an application to the observations requires a calculation of T3 for a field which has been 
smoothed with a window function. Both Fry & Scherrer (1994) and Chodorowski & Bouchet 
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(1996) argued that the smoothed skewness and kurtosis for non-Gaussian initial conditions should 
qualitatively resemble the unsmoothed results. The effects of smoothing on the hierarchical 
amplitudes for Gaussian initial conditions have been calculated in detail for the case of spherical 
tophat smoothing (Bernardeau 1994a,b), so we will follow Bernardeau's treatment to derive an 
expression for the value of T3 with non-Gaussian initial conditions and spherical tophat smoothing. 
In addition, we extend the results of Fry & Scherrer (1994) by performing the same calculation to 
derive the smoothed value of 53. 

Consider a spherical tophat window function with radius Rq and Fourier transform 

W ( kR o) = 77-5-^ [sm(kRo) ~ kR cos(kR )] . (35) 

After the density and velocity-divergence fields have been smoothed with this window function, 
we obtain new expressions for S3 and T3, which we denote by Ss(Ro) and T3(i?o). We now use the 
methods developed in Bernardeau (1994a) to calculate these quantities. The expressions for the 
smoothed 9 up to second order are (Bernardeau 1994a) 

j^S*W(kRo) 

9^(Ro) = -/(») / (2 ^ /2 { ^ /2 S kl S k2 W(\^ +k 2 |^) 

where <5k is the Fourier transform of the initial density field, the quantities Pij, Qij are defined by 

k? * k^' 



Dl(P 12 - \Qi2) + \E2Q12 



(36) 



Pij = 1 + 



k? 



Qij — 1 1.2 1,2^ ' ( 37 ) 

i 

and Di gives the time dependence of the growing mode, being the solution of the ith order time 

evolution equation for 5, while E 2 = Di——(D 2 — D 2 ). At the limit of t — ► 0, D 2 satisfies the 

dD\ 

34 

relation D 2 (t) « — D 2 (t). Furthermore, the 2,3 and 4-point functions of the random field <5 k are 
expressed as 

(£ki<*k 2 ) = £k 12 , (38) 
(<5 kl <5 k2 4 3 ) = Ck 123 , (39) 

((5 kl (5 k2 (5 k3 5 k4 ) = £ kl2 £ k34 + Ckl 3 Ck 2 4 + £k 23 £ki4 + »7kl 23 4> ( 40 ) 



where the 2-point function is related to the power spectrum P(k) of the (5 k field as 

£ k .. = 5 D (ki + kj)P(ki). Using these relations and given that (9 2 ) = ^ + 2 (6^6^) + 

we obtain 

[e 2 (Ro)) = f(Q) 2 D 2 a 2 (R ) + 2f(n) 2 Dfh [C kl23 ] (41) 
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where we have used the definitions 

2 



d 3 k 
(2vr)3 



W(kR ) 2 P(k), 



and 



-^1 [Cki 23 ] 



d 3 k ± d 3 k 2 d 3 k 3 

(2vr) 3 /2 ( 27r )3/2 ( 2yr )3/2 



W 1 W 23 ( kl23 



(P23 - 2^23) + 4 £j|<323 



where Wj = W(kiRo), and Wjj = W(|kj + \tj\Ro). In a similar fashion, we obtain 



0(1) 



= -mfD\i 2 [c kl23 ] , 



where we define 



h [Ckl2s] = 7 (2vr)3/2 ( 27 r)3/2 (2vr )3/2 W i^W3Ck 12 3. 



Finally, as shown in the Appendix, we get 



awl v> 



where we define 
Since (6 3 ) = 



d 3 h d 3 k 2 d 3 k 3 d 3 k 4 
(2vr) 3 / 2 (2vr) 3 / 2 (2tt) 3 / 2 (2tt) 3 / 2 



w x w 2 w 3m ^ 



0(1) 



+ 3 



0(1)1 0(2) \ + we obtain 



26 



(^34 - 2<334) + ^^1^34 



[e 6 (Ro)) = -fm 6 D 3 i 2 [c kl23 ] - /(^z^iio) (^y + 7 J - may mi [% 1234 ] 

where 7 is defined as 



7 



dlog [a 2 (R ) 



dlog(i? ) 



Using these results it is easy to show then that 

1 h [Ck 123 ] 



T 3 (R ) = - 



f(Q) D 1 a^(R ) 

3 / [%1234] , 4 A [Ck 123 ]^2 [Ck 123 ] 



+ 



1 f26 



f(Q) a\R ) f(Q) a^Ro) /(«) V 7 



(42) 
(43) 

(44) 
(45) 

(46) 
(47) 



(48) 
(49) 



(50) 
(51) 



As in the unsmoothed case, the first term decays as 1/D\{t), while the rest of the expression is 
constant in time. Observe that for Gaussian initial conditions our result reduces to 



1 /26 



m) v 7 



+7 , 



"26 



ft - 6 



- (n + 3) 



(52) 
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in agreement with Bernardeau (1994a), where the last relation holds for a power spectrum 
P(k) oc k n . 

Using the same line of argument, it is easy to show that for the smoothed density field, the 
mean values of 5 2 and 5 3 are 

;^-(i2o))=D?(7 2 (12 ) + 2^[^[u i ,:,]+ fv[c, kl2:) ]]. (.-,.■!) 
and 

(5 3 (Ro)) = D\I 2 [C kl23 ] + Dja\Ro) (y + t) + 3Df [J [r? kl234 ] + if [^ 1M4 ]] . (54) 

Then the skewness of the density field smoothed with a spherical tophat window function for 
non-Gaussian initial conditions is 



S 3 (Ro) 



h [Cki 23 ] 



^i^ 4 (^o) 
/34 \ 3 
+ (^y + 7j + ^7^j t J [^1234] + # [%i 23 J] 



[A [Ck 123 ] ^2 [Ck 123 ] + J 2 [Ck 123 ] # [Ck 123 ]] , (55) 



where we use the definitions 



jy- i 2 f f f d3k i ^2 d3k 3 n w w ^ 

K[CklJ = 7j J J (2vr)3/2 (2vr)3/2 Q^W 2 ^ 123 



(2vr) 3 / 2 (2vr) 3 / 2 (2vr) 3 / 2 
irr 1 2 f f f f rf3fcl d3fc2 d3fc 3 d3fc4 n WWW ^ 

^k 1234 ] = - 7 j j j j (27r)3/2 (27r)3/2 (27r)3 /2 ^3/2^1^34^ (56) 

and we have kept terms up to 0(<r°). For Gaussian initial conditions our result reduces to 

34 

S^(Ro) = y+7, 

= y-(n + 3), (57) 
for a power spectrum P(/c) oc k n , consistent with Bernardeau (1994a). 



4. Discussion 



A calculation of this sort can be used in two different ways: to gain further insight into the 
process of gravitational clustering at a fundamental level, and to compare with observations to 
try to determine if the initial perturbations were non-Gaussian. For the latter purpose, only our 
smoothed results are useful, although the unsmoothed results may be of interest for the former. 
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As noted earlier, the hierarchical amplitude for the velocity divergence, T3 resembles 
qualitatively the corresponding expression for the density field, S3. We end up with one term, 
proportional to the initial skewness, which decays as 1/Di, while the term which is constant in 
time can be broken down into a "Gaussian" piece, equal to the contribution for Gaussian initial 
conditions, and a "non-Gaussian" piece, which depends in a complex manner on the three- and 
four-point functions of the initial density field. 

A problem with applying the results of either Section 2 or Section 3 is their complexity 
in comparison with their Gaussian counterparts. Even the unsmoothed results are non-local, 
depending on integrals over the initial distribution functions. However, it is possible to simplify 
these results somewhat by examining combinations of S3 and T3. Consider first the unsmoothed 



case for = 1. In this case, the non-local terms arise from the last term in equation (12). These 
can be eliminated by evaluating 53 + 15T3, for which we obtain 

1 1 Co(0) wo) c 2 (o) 

s " + 2 T3 -2D 1 m + + W)~wr m 

This expression is a function only of the initial skewness and kurtosis of the non-Gaussian density 
field. 

A more useful combination, from the point of view of the observations, can be derived from 
the smoothed hierarchical amplitudes given in the previous section. For the case = 1, if we 
simply take the sum of Ss(Ro) and T^(Rq), we obtain: 

<7 (Ft \±T (R \ - 8 JU ^JWtmJ 4/ 2 [Ck 123 ]^ [Ck 123 ] r , Q x 

Ss{Ro) + Ts(jRo) - 7 + ^) — • (59) 

This quantity has several interesting properties. For the case of Gaussian initial conditions, it 
reduces to 

S 3 3 (R ) + T 3 3 (R ) = -, (60) 

which is independent of the initial power spectrum. This result does not extend to higher-order 
amplitudes; e.g., S4 — T4 does depend on the initial power spectrum. For non-Gaussian initial 
conditions, the time-dependent term produced by the initial skewness has vanished, so equation 
( |59| ) gives a much cleaner estimate of the deviation from Gaussian hierarchical clustering for 
non-Gaussian initial conditions; any deviation from 8/7 indicates the presence of non-Gaussian 
initial conditions or / 1. 

One of the main reasons for investigating the behavior of Ts(i?o) is its sensitivity to different 
values of Q (Bernardeau et al. 1995; Bernardeau et al. 1996). Unfortunately, this works to our 
disadvantage in equation (|||): it is not possible to disentangle the effects of 7^ 1 from the effects 
of non-Gaussian initial conditions. A more useful quantity if one is interested in the statistics 
of the initial conditions is ss(Ro) + t^(Ro), since £3 is independent of /(fi). A straightforward 
calculation similar to our derivation of Ss(Ro) + T^Rq) leads to: 

% K [%1234 ] 3/ 2 [Ck 123 ] K [Ck 123 ] 



S3(Ro)+t3(R ) 



+ 



7 a\R Q ) a%R ) 



DMR ), (61) 



a result which holds for any value of S7. Thus, equation (|6l"|) provides a clear distinction between 
the evolution of Gaussian and non-Gaussian initial conditions. 



The application of hierarchical amplitudes to the case of non-Gaussian initial conditions 
continues to be difficult, due primarily to the much greater complexity of the results. However, 
our calculations provide some simpler expressions which represent a step in the right direction. 
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Appendix 



In Section 2 we used the integrals 

1 



(47T) 



(4*) 



' //A<iV«*<-x<<)( * ),,( t^t, h 



= {(0), 



which are easily derived using an integration by parts and the relation 

1 3xj • Xj - Sij 4tt 



x 



(62) 
(63) 

(64) 



In Section 3 we used the results given by Bernardeau's equations (A26) and (A27) (Bernardeau 
1994a), specifically, 



d 3 h d 3 k 2 



[ [ d k 
J J (2^ 



(2tt) 3 (2tt) 3 

d 3 h d 3 k 2 
(2tt) 3 (2nf 



P 1 2W 1 W 2 W 12 P(k 1 )P(k2) = a\R ) 1 + 
Q 12 W 1 W 2 W 12 P(k 1 )P(k2) = V(i?o), 



1 Rp da 2 (R ) 
6a 2 (Ro) dR 



(65) 

(66) 
(67) 



and the fact that in the equation for ^ 8^ 9^ y, integrals containing the factors Pjj<5(kj + k,-) 
and Qij5(ki + k,) vanish. 



